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Abstract 

We study the effect of going beyond the Gross-Pitaevskii theory on the frequencies of collec- 
tive oscillations of a trapped Bose gas in the large gas parameter regime. We go beyond the 
Gross-Pitaevskii regime by including a higher-order term in the interatomic correlation energy. To 
calculate the frequencies we employ the sum-rule approach of many-body response theory coupled 
with a variational method for the determination of ground-state properties. We show that go- 
ing beyond the Gross-Pitaevskii approximation introduces significant corrections to the collective 
frequencies of the compressional mode. 

PACS numbers: 03.75.Fi, 03.65.Db, 05.30.Jp 
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I. INTRODUCTION 



The mean-field Gross-Pitaevskii (GP) theory Jl| has been quite successful in explaining 
both static and dynamic properties [Q, 0, |4|] of Bose-Einstein condensates (BEC) produced in 
alkali atoms [[5], H [7|] • The reason for the success of the GP theory is mainly the satisfaction 
of dilute gas condition na 3 << 1 ( where n is atomic density, a is the s-wave scattering 
length of interatomic potential and the parameter na 3 is called the gas parameter) in the 
above mentioned experiments. Typically the values of gas parameter in these experiments 
were in the range of 10 -5 — 10~ 4 . The gas parameter can have a large value when the 
number of atoms N in the condensate or the scattering length is large. Thus in the later 
experiments reported in Refs. || ||| N was of the order of 10 s and the corresponding 
value of gas parameter was 10~ 3 . On the other hand, in a very recent experiment |10] the 
maximum value of the peak gas parameter was n(0)a 3 ~ 10~ 2 (^(0) is the peak density of 
the condensate) although the condensate contained only 10 4 85 Rb atoms. In this experiment 
the scattering length a could be varied continuously from a negative to a very large positive 
value by exploiting the Fesbach resonance. These kind of condensates, then, naturally raise 
the question about the validity of the mean-field GP approach. Theoretical investigations of 
the effects of going beyond the GP (BGP) theory have already been reported in the literature 
fTT| , |T!| |T3| , |T5l . Recently []16[ ground-state properties of BEC with large values of gas 
parameter (as achieved in the experiments (1(|) have been theoretically studied to assess 
the accuracy of the GP theory. It has been shown that for large values of the gas parameter 
(« 10 -2 ) the total energy and the chemical potential are modified significantly by the BGP 
theory. It is well known that these changes are also reflected in the frequencies of collective 
excitations which can be measured with a greater accuracy. Therefore, by measuring the 
frequencies of collective excitations it should be possible to observe the effect of the BGP 
theory. This has motivated us to calculate the corrections to the frequencies of collective 
oscillations due to the BGP theory. We focus our attention on the condensates with values 
of gas parameters, spanning a range similar to that achieved in the experiment [|TUJ and make 
an assessment of the accuracy of the GP theory in predicting the frequencies of collective 
oscillations of such condensates. 

In the past corrections to the frequencies of collective oscillations of trapped bosons aris- 



ing due to the BGP theory have been estimated JT7|, [TJ| within the Thomas-Fermi (TF) 



2 



approximation which is valid in the large N limit (more precisely when Na/ciho » 1 is sat- 
isfied, where a^o is the characteristic length of trapping potential and it is defined in the next 
section). In the aforesaid references corrections to the frequencies were obtained by solving 
the time-dependent second-order differential equation perturbatively for the density. Thus 
the analytical expressions for the fractional change in the frequencies are correct only up to 
the first order in ^n(0)a 3 . It is, then, expected that these expressions for the corrections 
will not be accurate in the large-gas-parameter regime, as achieved in the recent experiment 
Til . In this paper we also derive an analytical expression for the correction within the TF 



approximation by making use of the sum-rule approach. We find that the expression for the 
correction obtained by us is valid for somewhat larger range of gas parameters than those 
of Refs. |L7], [18| . As expected, for very small gas parameters the expression derived by us 



correctly reduces to those of Refs. |17], |18| . We will return to the question of validity of 
these analytical results later. 

It is well known that the TF approximation does not reproduce the surface region of 
density profile of trapped condensates accurately [Q]. On the other hand, the collective 
oscillation frequencies characterizing the response properties of the trapped gas crucially 
depend on the tail or the surface region of the density profile |19[]. Therefore, it becomes 
necessary to go beyond the TF approximation especially for the calculation of frequencies of 
collective oscillations [13]. To accomplish this task we make use of well established variational 
method [T^, |2D|, to determine the ground-state densities. The collective oscillation 

frequencies are then obtained by using the sum-rule approach |23|, |24| of many-body response 
theory. In the sum-rule approach the calculation of frequencies requires the knowledge of 
the ground-state wave function (or the ground-state density) of the many-body system only. 
Therefore the accuracy of frequencies determined by the sum-rule approach depends on 
the accuracy of the ground-state wave function (ground-state density) employed for the 
calculations. Keeping this in mind our aim is also to obtain reasonably accurate ground- 
state densities of the trapped condensates. To this end we employ variational approach 
by using an energy functional of N interacting trapped bosons within the local density 
approximation (LDA) along with a judiciously chosen variational ansatz for the ground- 
state density [15|, |22| . 

The paper is organized in the following manner. In Section II we describe the theoretical 
method for the calculation of the frequencies of collective oscillations. Section III contains 
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the analytical estimate of the corrections to the frequencies within the TF approximation, 
and results of our variational calculation. The paper is concluded in Section IV. 



II. VARIATIONAL SUM-RULE APPROACH 

The central quantity in the variational approach is the energy functional for a condensate 
of iV bosons, each with mass m, confined in a trap potential Vj(r): 



E\n] = / dr 
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- — A/"-( r )V 2 A/n(r) + V t (r)n(r) + e(n)n(r) 



(1) 



where, the first, the second, and the third terms represent the kinetic energy of bosons, the 
energy due to the trapping potential, and the energy due to interatomic correlation within 
the LDA, respectively. Following the usual practice we use the perturbative expansion for 
e(n) in terms of the gas parameter na 3 , 



e(n) 



2nh an 



1 + yi^= (na 3 ) * + 8 - V3j (na 3 ) In (na 3 ) + O (na 3 ) 



(2) 



m 

The first term in the above expansion, which corresponds to the energy of the homogenous 
Bose gas within the mean-field theory as considered in the GP theory, was calculated by 



Bogoliubov [25|. The second term was obtained by Lee, Huang and Yang (LHY) [pq] , while 
the third term was first calculated by Wu |27] using the hard-sphere model for the interatomic 
potential. Although, it has been emphasized in the literature that the above expansion is 
valid only for na 3 << 1, it is only recently that the range of validity of above expansion has 
been systematically investigated by Giorgini et al. ||28|| . They have used diffusion Monte 
Carlo (DMC) method to calculate the ground-state of uniform gas of bosons interacting 
through different model potentials with a > 0. It has been found that the expansion (^|) is 
valid as long as na 3 < 10~ 3 . However, for the values of the gas parameter beyond 10 -3 the 
inclusion of the logarithmic term leads to severe mismatch with the DMC simulation results. 
On the other hand, expansion @ up to the LHY term gives an accurate representation of 
the DMC calculations even for the gas parameter of the order of 10~ 2 . Consequently, we do 
not consider the logarithmic term in the expansion (0) for all the calculations in this study. 

The trapping potential V t (r) is taken to be axially symmetric characterized by two angular 
frequencies u Q ± and u° z (cu° = = ^ u®). It is given by 

V t (r) = ^f(xZ + y* + \y), (3) 
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where Ao = uj^/uj ^ is the anisotropy parameter of the trapping potential (Aq = 1 corresponds 
to a spherically symmetric trap). 

The ground-state properties which are needed to calculate the collective frequencies are 
determined by minimizing the energy functional ([I]) with respect to the density n(r) 
satisfying the normalization condition 



J n{r)dr = N. (4) 



To implement the variational scheme we take the ansatz for the ground-state density as 



3 / \ P 

2 I \ (r?.-l-\x?Y 



n(riH ^(4) A2 UJ eUy ' (5) 

where A, u± and p are the variational parameters which are obtained by minimizing E[n] 

with respect to these parameters. In the above equation we have used the scaled co-ordinates 

i 

i"i = r/dho with a/j = (h/mu°^) 2 . We note that n(ri) is normalized to unity and from Eq. 
(|) it can be seen that 

N 

n(r) = — n(ri). (6) 

a h0 

We make this choice of variational form for the density as it has been shown convincingly that 
this ansatz is capable of describing the ground state of Bose gas in a trap quite accurately 
both within |22| and beyond [[HJ the GP theory for a wide range of particle numbers. Using 



this variational form the physical observables can be expressed analytically in terms of three 
variational parameters ]T5|] . The total energy Ei = E/fiWj_ can be written as 

^ = T + U + El nt + El t . (7) 

Here T and U denote the average kinetic and the trapping potential energies per particle, 
respectively. The term E] nt gives the interaction energy per particle in the mean-field ap- 
proximation as considered in the GP theory while Ef nt gives the correction in the interaction 
energy due to the LHY term in the expansion (H) arising from the BGP theory. 

The ground-state energy components obtained from the variational calculation are em- 
ployed to compute the frequencies of collective oscillations by using the sum-rule approach 
of many-body response theory. In the following we briefly outline the sum-rule approach 
and present some of the results which are relevant for the present work. 



According to the basic result of the sum-rule approach |23|, ^4] the upper bound of the 
lowest excitation energy is given by 



hQ ex = J^, (8) 
V mi 

where 

^\(0\F\n)\ 2 (hw nQ r, (9) 



m p 



is the p-th order moment of the excitation energy associated with the excitation operator 
F and Q ex is the frequency of excitation. Here hu n0 = E n — E is the excitation energy of 
eigenstate |n) of the Hamiltonian H of the system. The upper bound given by Eq. (^) is 
close to the exact lowest excited state when this state is highly collective, that is, when the 
oscillator strength is almost exhausted by a single mode. This condition is satisfied by the 
trapped bosons in most of the cases. Moreover, Eq. ([8]) can be used for computation of the 
excitation energies by exploiting the fact that the moments mi and can be expressed as 
expectation values of the commutators between F and H in the ground state |0) |^3. [2~4|j : 



mi = ~<0|[Ft,[F,F]] |0>, 

m 3 = ~(0\[[F\H],[[H,[H,F]}]\0). (10) 

The main advantage of the sum-rule approach is that it allows us to calculate the dynamic 
properties like excitation frequencies of many-body systems with the knowledge of ground 
state |0) (or the ground-state density) only. Thus, with the knowledge of a reasonably accu- 
rate ground-state density we calculate the expectation values given in Eq.flTOD to determine 
the frequencies of collective oscillations. It is worth mentioning that sum-rule approach has 
been extensively used to calculate collective frequencies of the trapped condensate pHJ. 

For the purpose of calculation one needs to choose an appropriate excitation operator 
F. In this paper we concentrate only on the compressional collective mode characterized by 
the z-component of angular momentum index, m' = 0. For the axially symmetric trap the 
mode ml = involves coupling of the monopole (I = 0) and the quadrupole (/ = 2) modes. 
Following Kimura et al. |3l| the operator F for the excitation of m' = collective mode is 
written as 

F = Y,{^ + yl-azl) ) (li) 
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where a is a parameter characterizing the coupling of the two modes due to the axial 
symmetry of the trap. This parameter is obtained by making the excitation energy given by 
Eq.(^|) extremal. Here, we note that F given above can also be used to describe excitation 
of the monopole (a = —1) and the quadrupole (a = 2) modes separately. For spherically 
symmetric trap (Ao = 1) the two modes get decoupled and a = — 1 and a = 2 correspond 
to the minimum and the maximum values of the excitation energies, respectively. 

Using the energy functional given by Eq. (jH) along with the expansion @ we find after a 
tedious although straightforward algebra the following expressions for the moments mi and 
m 3 : 

Ah 2 I r 



mi 



10J 



0\* 



+ a" 



\2 I ' 



(12) 



m 3 



8h 4 



(13) 



In the above equations T± and ?7j_ ( T z and ) are the transverse components (z- component) 
of the kinetic and potential energies, respectively. By substituting Eqs. (0) and ( p~3|) in Eq. 
(§) we get following expression for the frequency of m! — mode 

_ JT± + U ± + o? (T, + U z ) + (1 - a/2) 2 (E} nt + \E^)] 



m'=0 



C7 



_|_ n-2 ^ 



(14) 



The frequency of collective oscillations is determined by substituting the values of energy 
components from the variational calculation and making fi^/ =0 extremal with respect to 
a. The variation with respect to a leads to a quadratic equation in a and the two roots, 
corresponding to the maximum (a+) and the minimum (a_) values of f2^, =0 , are 



with 



B 



(2/2/4 + J3/4/2 - 2/i/ 5 - 2/ 3 / 5 



(15) 



and 



C 



h 



(16) 
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where fi = T ± + U ± , f 2 = T z + U z , f 3 = E} nt + \E\ nt , / 4 = U ± /(u° ± ) 2 and f 5 = U z /(u° z ) 2 . 
In the next Section we first apply these results to estimate analytically the correction to 
the collective frequency due to the LHY term within the TF approximation followed by the 
results of our variational calculation. 



III. RESULTS AND DISCUSSION 



A. Analytical estimate within TF approximation 



Before discussing the results of our variational calculation we present the analytical esti- 
mate of collective frequencies within the TF approximation. We generalize the perturbative 
results of Ref. |T7|, [18] and also discuss the validity regime of these results. For convenience 
we consider the case of spherically symmetric trap (Ao = 1). Generalization for the axi- 
ally symmetric trap is straightforward. It is easy to verify that for spherically symmetric 
trap B = — 1 and C = 2 and the two roots are a = 2 and a = — 1 corresponding to 
the quadrupole and the monopole modes, respectively. It is well known that the frequency 
of quadrupole mode is not affected by the LHY term |17| |. Therefore we focus our atten- 
tion on the monopole mode. By substituting a = — 1 in Eq. fll4]) the frequency of the 
monopole mode of bosons in a spherical trap (U = U± + U z ) within the TF approximation 
(Tj_ = T z = 0) can be written as 



Q 2 



\ 2 



u 



(17) 



Hereafter we drop for convenience the subscript m' = in Q. To eliminate U from the above 
equation we make use of the virial relation 



2T - 2U + 3El 



int 



-E 



int 



0. 



(18) 



and obtain within the TF approximation 



Q 2 



[UJ, 



15 E l 
4 int 



16 int 



3 rpl , 9 rp2 
2 JZj int 4 int 



(19) 



We note that the expression given by Eq. (|T9|) is exact within the TF limit. To derive 
an analytical expression for the frequency we express E} nt and E 2 nt in terms of parameter 
x = y / n(0)a 3 . This is done by using the analytical expressions for the total energy and 



S 



chemical potential 0, |15| obtained within the TF approximation with the LHY term included 
in the expansion e(n). These expressions are: 

pi 2 9 

EL = I^tfx, (20) 

where \xtf is the chemical potential of trapped bosons obtained within the TF approximation 
0. It is important to note that the expressions for E\ nt and Ef nt given above (Eq. (|2"0D) are 
approximate in nature and neglect terms of order 0(na 3 ) and, thus, expected to be valid in 



the small gas parameter regime. It has been recently shown |16| that the expressions in Eq. 
( p0|) are not valid for very large gas parameters and lead to results which differ considerably 
from the exact results . In the large-gas-parameter regime one needs to solve the nonlinear 
equation of the density more accurately with the help of numerical methods. Consequently, 
to compute the frequencies of collective oscillation (for that matter any observable) in the 
large-gas-parameter regime within the TF approximation one needs a non-iterative or non- 
perturbative method to solve the nonlinear equation. In this paper we take a different route 
to achieve this goal. Instead of solving the time-dependent differential equation for the 
density (or the time-dependent GP equation) we exploit the sum-rule approach of the many- 
body response theory along with the accurate ground-state density determined variationally 
to compute the frequencies of collective oscillations. We describe the results of such a 
calculation in the next subsection. 

Now by substituting Eq. fl20|) in Eq. ([191) we get following expression for the frequency 
of the monopole collective mode within the BGP theory: 

1/2 

(21) 



LO 



5 (1 + W*) 



where uq is the angular frequency of the spherically symmetric trap. In the limit of x — > 
above expression (Eq.(|2~lf)) correctly reduces to the mean- field result VL m f = \Z5to within 
the TF approximation (3( 



It is customary to define the fractional change in the frequency as |L7], [18 
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By using Eq. fl21|) we get 

63 



l2s '' (23) 



In the limit of x << 1 we recover the results of Refs. |T7], |Tj|. Fig. 1 shows the 
fractional change 6 as a function of dimensionless parameter a/ao (where ao is the Bohr 
radius of hydrogen atom) for iV = 10 4 85 Rb atoms trapped in a spherical trap with 
uj /2it = (V5J 2 wf) 1/3 /27r = 12.83Hz, where u° ± = 2vr x 17.5Hz and lu° z = 2vr x 6.9Hz 
fT0|| . The values of the parameter a/ao are chosen such that the condition for the TF 
approximation, that is , Na/a^Q » 1 is satisfied. Moreover, as discussed above Eq.(pO|) 
(and hence Eq. fl2"T|)) is valid only for small gas parameter values and accordingly we have 
restricted ourselves to a/ao = 5000. It is evident from Fig.l that the correction to the 
monopole frequency due to inclusion of LHY term in the interatomic correlation energy 
obtained by us is slightly lower than the corresponding results of Refs. ||17| , |18|| . In partic- 
ular for a/ao = 5000 we obtain 7% change in the monopole frequency as compared to the 



correction of 9% predicted in Ref. [|T7| , [18]. It is only for very small values of x (less than 
10 -3 ) that fractional corrections obtained in Ref. [17, 18] match with Eq. (|23|). 



B. Variational results 

Before discussing the results for the excitation frequencies obtained by variational method, 
we first demonstrate that this method is capable of yielding sufficiently accurate results for 
the ground-state properties of BEC even for large gas parameter values, as achieved in the 
experiment of Cornish et al. ||10|| . The parameters used in this calculation are similar to 



the ones used in the previous section. We present these results in Table I and II for the 
spherically symmetric and the axially symmetric traps respectively. To check the accuracy, 
we compare our variational results with those of Ref. |16[ (presented in parenthesis in Table 



I and II) obtained by numerically solving the modified GP equation. The results in Table I 
and II clearly show that the variational method gives quite accurate results for the ground- 



state properties as they compare very well with the numbers of Ref. ||16|| . Therefore the 
variational method is indeed suitable for describing the ground-state properties even in the 
regime of large gas parameters. The accuracy of our variational results was further checked 



by verifying that the generalized virial relation as given in Eq. (|TB"D is satisfied. In all 
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the variational calculations the virial relation was satisfied up to the fifth decimal place or 
better. This indicates that not only the total energy but its different components are also 
obtained accurately by the variational method, which is quite crucial for the calculation of 
the frequencies of collective oscillations. 

Now we focus our attention on the results for the frequencies of collective oscillations 
of trapped bosons. Figs. 2 and 3, show the upper (Q u ) and the lower (fij) frequencies 
respectively, associated with the breathing mode as a function of the dimensionless parameter 
a/cio- The results presented in Figs. 2 and 3 correspond to the axially symmetric trap which 
is more relevant from the experimental point of view. Notice that the maximum value of the 
parameter a/a^ is chosen to be 10 4 in accordance with the experimental number [ 10 1 and 
the corresponding value of the peak gas parameter is ~ 10 -2 . For comparison we also show 
in these figures the corresponding results obtained within the GP theory. It is important 
to note here that the GP results shown in Figs. 2 and 3 are also obtained variationally by 
dropping the LHY term from the expansion (0) and thus, retaining only the first term. The 
results corresponding to the GP case obtained by us match quite well with those already 



reported in the literature [30, 31, 32, 33, 34, 35, 36, pTl 38|. It can be clearly seen from 



Fig. 2 that for small values of the gas parameter the numbers obtained within the GP 
and the BGP theory are quite close. The difference between these two results, however, 
increases with a. For example, for a/a = 8000 the BGP number for Q u is 6% higher than 
the corresponding mean-field GP result. Thus we conclude that the introduction of LHY 
term in the interatomic correlation energy results in a significant correction to the upper 
frequency Q u of m! — mode of the collective oscillations. Considering the accuracy of 
current experiments |§ in measuring the frequencies of the collective excitations, it should 
be possible to observe these corrections experimentally. It is also important to note that 
the frequency Q u obtained within the GP theory remains almost constant with the increase 
in the parameter a/oo- In contrast, the BGP theory predicts Q u growing monotonically 
with a/ do- Therefore a signature of the predictions of the BGP theory can be observed by 
measuring a/a^ dependence of the excitation frequency fl u . 

Next we focus our attention on the lower frequency fli of the breathing mode which is 
shown in Fig. 3. For small values of the gas parameter the frequency of the lower branch 
obtained via GP theory is very close to the corresponding BGP theory result. In the same 
regime of the gas parameter the frequency Qi decreases with the gas parameter both for the 
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GP and the BGP cases. The decrease of the lower frequency with an increase in the gas 
parameter value has already been theoretically predicted 0, |32], |57] and experimentally 
verified as well |]3*9fl . However, in contrast to the GP case, frequency of the lower branch 
for the BGP case starts increasing, albeit slowly, after a critical value of a/ao ~ 2200. 
It is also noteworthy that the change in the frequency Qi obtained from the BGP theory 
is considerably lower than that in the upper frequency Q u . For example, for N = 10 4 and 
a/ao — 8000 the change in the lower frequency is just 0.8% while that in the upper frequency 
is 6%. 

We conclude this section by making a comparison of the frequencies obtained by using 
the analytical expression as give in Eq. (|2T|) with the results obtained by the variational 
approach. The comparison is shown in Fig. 4 where we plot the frequency of the monopole 
mode of a condensate in a spherically symmetric trap as a function of a/ao obtained using 
Eq. (|21~D (dashed line) and using the variational method (solid line) . The range of values of 
a/ao is chosen such that the TF approximation (i.e. Na/aho >> 1 ) is satisfied. Accordingly 
the minimum value of a/ao m Fig. 4 is chosen to be 10 3 which corresponds to Na/aho — 174. 
It can be seen that the results obtained using Eq.(^l|) are systematically higher than the 
corresponding results from variational calculation and the difference between the two grows 
as the value the gas parameter is increased. Moreover, we also note that the two results are 
reasonably close for values of a/ao of the order of 3500 (ra(0)a 3 ~ 3.77 x 10~ 3 ). Therefore, 
we emphasize that the correction obtained within the TF approximation as given by Eq. 
(^T|) is valid as long as the gas parameter is of the order of 10~ 3 . As discussed earlier, for 
condensates with the values of the gas parameter beyond 10 -3 one needs to evaluate E\ nt 
and Ef nt in a non-perturbative manner. The variational approach employed in this paper 
accomplishes exactly this task. 

IV. CONCLUSION 

In this paper we have calculated the frequencies of collective excitation of a trapped 
condensates in the large-gas-parameter regime using the BGP theory . To go beyond the 
GP theory we have included a higher-order term (LHY term) in the interatomic correlation 
energy obtained from the ground-state energy of a uniform Bose gas. The calculations of 
the frequencies were performed by employing the sum-rule approach of many-body response 
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theory along with the variational method to obtain the ground-state energy components. 
We have also derived an analytical expression for the corrections to the frequencies of the 
compressional mode valid for small gas parameters within the TF approximation. The 
results of our variational calculations clearly show that for large gas parameter values, the 
frequencies of the compressional mode characterized by the index m' = are significantly 
altered. Although the change in the lower branch of m' = mode, f^, is small, we observe 
a new qualitative feature. The frequency Qi increases slowly as the gas parameter exceeds 
a critical value in contrast to a monotonically decrease predicted within the GP theory. It 
would be very interesting to verify the predictions of the BGP theory by measuring collective 



excitation frequencies for condensates with large gas parameter as reported in Ref. [10 
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TABLE I: Results for the ground-state properties of Rb atoms trapped in an isotropic trap with 
lJ\_/2it = Yl.hHz. Total energy E\/N and chemical potential fii are in units of fko^_. Numbers in 
parenthesis are results of Ref. 113]. The values of parameter Na/a^o given in the second column 
signify the importance of the interaction energy over the kinetic energy. 







Ei IN 


pi 


i4nn 


243 38 


1 03 


14 05 






(10.01) 


(13.97) 


3000 


522.0 


14.14 


19.92 






(14.10) 


(19.84) 


8000 


1392.0 


23.52 


33.54 






(23.43) 


(33.38) 


10000 


1740 


26.75 


38.24 






(26.63) 


(38.06) 



TABLE II: Results for the ground-state properties of 85 Rb atoms trapped in axially symmetric 
trap with uP^jl-n = 17.5Hz and Ao = 0.39. Total energy E±/N and chemical potential [i\ are in 
units of hu>^_. Numbers in parenthesis are results of Ref. J16| |. The values of parameter Na/a^o 
given in the second column signify the importance of the interaction energy over the kinetic energy. 



a/a 


Na/a h0 


Et/N 


Ml 


1400 


284.25 


7.36 


10.25 






(7.33) 


(10.22) 


3000 


609.0 


10.35 


14.55 






(10.31) 


(14.51) 


8000 


1624.0 


17.18 


24.49 






(17.09) 


(24.38) 


10000 


2030.0 


19.53 


27.92 






(19.43) 


(27.79) 
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FIG. 1: Fractional shift of the frequency of m' = mode in a spherically symmetric trap with 
loq = 2tt x 12.83-ffz. The solid line corresponds to Eq. ( |23| ) and the dashed line shows the results 
of Fiefs. 0,0. 



FIG. 2: The frequencies (in units of uj^ ) of the upper branch of the collective mode m' = of 10 4 
85 Rb atoms confined in an axially symmetric trap with Ao = 0.39 and = 2ir x ll.hHz. The 
solid line represents results for the BGP case while the corresponding GP results are shown by the 
dashed line. 

FIG. 3: The frequencies (in units of lu^ ) of the lower branch of collective mode rn! = 0. The 
parameters are same as in Fig. 2 

FIG. 4: The frequencies (in units of luq ) of the monopole mode of 10 4 85 Rb atoms confined in a 
spherically symmetric trap with parameters same as Fig. 1. The solid line represents results for 



the BGP case and the dashed line shows results obtained from Eq. (21) 
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